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Abstract. We show that the amount of coherent quantum information that can be 
reliably transmitted down a dephasing channel with memory is maximized by separable 
input states. In particular, we model the channel as a Markov chain or a multimode 
environment of oscillators. While in the first model the maximization is achieved for 
the maximally mixed input state, in the latter it is convenient to exploit the presence of 
a decoherence-protected subspace generated by memory effects. We explicitly compute 
the quantum channel capacity for the first model while numerical simulations suggest 
a lower bound for the latter. In both cases memory effects enhance the coherent 
information. We present results valid for arbitrary size of the input. 
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1. Introduction 

Quantum communication channels [U [2] use quantum systems to transfer classical 
or quantum information. In the first case, we can encode classical bits by means of 
quantum states. In the latter case, we may want to transfer an unknown quantum state 
betweeen differents units of a quantum system, for instance of a quantum computer, or to 
distribute entanglement between communicating parties. In both cases, the fundamental 
question is what is the maximum rate of classical or quantum information that can 
be faithfully transmitted. Classical and quantum capacities, defined as the maximum 
number of bits/qubits that can be reliably transmitted per channel use, provide the 
answer to this question. 

Quantum channels with memory are the natural theoretical framework for the study 
of any noisy quantum communication system where correlation times are longer than 
time between consecutive uses. This scenario applies to optical fibers which may show 
a birefringence fluctuating with characteristic time longer than the separation between 
successive light pulses [3] or to solid state implementations of quantum hardware, where 
memory effects due to low-frequency impurity noise [4] produce substantial dephasing 



Some theoretical result on quantum channels with memory has been already 
discussed for transmission of both classical and quantum information through a quantum 
channel. With regard to classical information transmission down a memory channel, it 
was pointed out that it can be enhanced by using entangled input states [HI [3, [8], 
and coding theorems have been recently proved for classes of memory quantum 
channels [9j [10] . Concerning quantum capacity, a lower bound has been found for some 
classes of channels with memory [TT] and subsequently specific model environments 
(structured in two parts, one responsible for memory effects and the other acting as a 
memoryless environment) have been studied [T2 l [T3l [Ti] . In particular, coding theorems 
for quantum capacity have been proved in [TJ] for the so-called forgetful channels, for 
which memory effects decay exponentially with time. 

The problem is formalized by considering the iV-uses Hilbert space Ti^ = 7-[® N and 
defining the system S, described by the reduced density matrix (RDM) p for iV uses. 
The input state is p = YsfLiPiPi, namely states chosen from the ensemble {pi, px}, 
with a priori probabilities {pi, ...,p K }, are sent down the channel. Due to the coupling 
to further uncontrollable degrees of freedom, the transmission of S may be noisy. The 
output is therefore described by a linear, completely positive, trace preserving (CPT) 
map £n(p), corresponding to iV-uses (the single use is defined in H and described by 
£ ). The map £n(p) can always be represented starting from an enlarged vector space 
including a suitably chosen environment E, initially in a pure state: Wq = |0) E (0| 



where U is a suitable unitary evolution of S+E referring to N uses. The conditional 
(depending on p) evolution of the environment can also be considered. It is described 



0. 




(1) 
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by the environment RDM and allows to define the conjugate CPT map, w = Tr s [U (p<S> 
w )W] =:S N (p). 

The quantum capacity Q refers to the coherent transmission of quantum information 
(measured in number of qubits), and it is related to the dimension of the largest subspace 
of Ti,N reliably transmitted down the channel, in the limit of large N. The value of Q 
can be computed, for memoryless channels, as [151 EEl EE3 EH US] 

Q = lim Q N = m&xI c (£ N , p), (2) 

N^oo iv P 

I c (£ N ,p) = S[£ N (p)]-S?(p). (3) 

Here S(p) = — Tr[plog 2 p] is the von Neumann entropy, S^(p) = S[£n(p)] is the entropy 
exchange [20]. The quantity I c (£n,p) is called coherent information [21] and must be 
maximized over all input states p. 

The limit N — > oo in (j2]) makes difficult the evaluation of Q. On the other hand this 
regularization is necessary, since in general J c is not subadditive. Indeed for entangled 
input states p [16J we may have I c (£]^,p) > J2k=i ^c(£> where p^ = Trs_( fc )(p) 
refers to the individual transmission of the k— th unit of information, therefore in general 
it cannot be excluded that Qn/N > Qi- The regularization is not necessary if the final 
state w of E can be reconstructed from the final state p' of the system. In this case, 
referred to as degradable channels [221 [231 [2H [25], it exists a CPT map T such that 
£ = T o £. It turns out [22J that for degradable channels the coherent information 
I c (£n,p) reduces to a suitable conditional entropy pQ, which is subadditive and concave 
in the input state p, and therefore the quantum capacity is given by the "single-letter" 
formula Q = Q\. 

In this work we focus on dephasing channels with memory. Dephasing channels 
are characterized by the property that when N qubits are sent through the channel, 
the states of a preferential orthonormal basis {\j) = \ji, ....,j N ), ji, ...,j N = 0,1} 
are transmitted without errors, implying a conservation law to hold [26]. Therefore, 
dephasing channels are noiseless from the viewpoint of the transmission of classical 
information, since the states of the preferential basis can be used for encoding classical 
information. Of course superpositions of basis states may decohere, thus corrupting the 
transmission of quantum information. Dephasing channels are relevant for systems in 
which relaxation is much slower than dephasing [271 H] • When memory effects are taken 
into account, we have that En ^ £® N , i.e. the channel does not act on each carrier 
independently. 

We show that the coherent information is maximized by input states separable and 
diagonal in the reference basis In particular, we calculate the coherent information 

for two models of dephasing channels. For a Markov chain we show that the coherent 
information is maximized by maximally mixed input states and compute Q. For an 
environment modeled by a bosonic bath, we propose a coding strategy based on the 
existence of a decoherence-protected subspace generated by memory effects and use 
numerical results to suggest a lower bound for Q. It turns out that in both cases 
memory effects increase the coherent information. 
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2. The dephasing channel and quantum capacity 

The unitary representation of the generalized dephasing channel [22] reads 

C/|j)|0 E ) = |i)|^), (4) 

where \<f>j) are environment states, in general non mutually orthogonal, describing the 
conditional evolution. The map Sn can be written in the Kraus representation [H [2] as 

p' = S N (p)=J2A aP Al (5) 

a 

where the system operators (A a )ji = (a E \(fij) Sji are diagonal in the reference basis (here 
{|ct E )} is an orthonormal basis for the environment). It is easily shown that this channel 
is degradable [22]. Indeed, for a generic input p = J2j.i Pji\j) equation (pEJ) yields 

w = £ N (p) = ^Pjj\^j){<f>j\. (6) 

3 

Since w only depends on the populations pjj which are conserved, we can write as well 
£ n = £n ° £, thus proving degradability. 

We now show that for a generalized dephasing channel the coherent information 
I c (£n,p) is maximized by input states diagonal in the reference basis. To this end we 
introduce 

p k = l{p k - 1 + ^p k ^), (k = l,...,N), (7) 

where p = p and the local operator £( fc) = l (1) <g>- ■ ■<S)t {k ~ 1) ®a { z k) <g>l (k+1) <g)- ■ -(g)l (iV) acts 
non-trivially only on the k— th qubit, by the Pauli operator which has eigenvectors 
\jk). We can easily see that pn is the diagonal part of p, by using the standard 
representation of the iV-qubit density matrix: 

P = X! Ch...i N CTil ® • • • ® cr£°, i k = 0, x, y, z, (8) 

{ik} 

where o"o = 1. We now study the action of the operators S^. First of all 
£jv(£( fc ' 'pS^) = Tj^£jq{p)Tj^ for any k and p, since E^* 1 commutes with the Kraus 
operators in ©. Also S[E{ k ^£ N {p)E^} = S[£ N {p)}, since the von Neumann entropy is 
invariant under unitary local transformations. Moreover £^{11^ p£^) = £jv(p), since 
the populations of S^pS^ are the same as for p. We can therefore conclude that 
I c (£n, £Wp£00) _ j c (£ Ni p^. This latter relation, together with the concavity of the 
coherent information for degradable channels (a direct consequence of the concavity of 
the conditional von Neumann entropy) implies that 

I c (£n, Pn) > L{.£n, Pn-i) > ■ > I c {£n, Po)- (9) 

Hence, diagonal input states maximize the coherent information. These states are 
separable, since they can be written in the form 

PN= E Q31-.-3NP?? ® ■ ■ ■ ® P? N \ (10) 
3i,-,3N 

with pfj = \j k )(jk\, (k = 1, ...,N), < q h ... jN < 1 and Y j: , x <i, ... JS = 1. 
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3. The memory dephasing channel 

3.1. Forgetful channels 

Interesting results on the quantum capacity of dephasing channels with memory can 
be obtained for forgetful channels, for which the memory dies out exponentially with 
time. Forgetfulness is defined in [14], according to a model in which the environment 
is structured in two parts: a memoryless one and one responsible for memory effects 
(see also [12] ). A key feature of forgetfulness is that it permits, with a negligible error, 
the mapping of the memory channel itself into a memoryless one. This may be clarified 
by referring to the double-blocking strategy p3]: we consider blocks of iV + L uses 
of the channel and do the actual coding and decoding for the first N uses, ignoring 
the remaining L idle uses. The resulting CPT map £n+l acts on density matrices p 
on TL® N . If we consider M uses of such blocks, the corresponding CPT map £m{n+l) 
can be approximated by the memoryless setting (8(n+l))^ M . This is possible because 
correlations among different blocks decay during the idle uses. This property can be 
expressed as follows[14j: 

\\£~m(n + l)( P ) - (£ {N+ L)f M (p)h <h(M- l)c" L , (11) 

for any input state p in '}i® MN ; where c > 1, || ■ ||i is the trace distance PQ, and h is some 
constant depending on the memory model (note that c and h are independent of the 
input state). This equation states that, even though the error commited by replacing 
the memory channel itself with the corresponding memoryless channel grows with the 
number M of blocks, it goes to zero expontially fast with the number L of idle uses in a 
single block. Equation (fTTl) permits the proof of coding theorems for forgetful quantum 
memory channels, by mapping them into the corresponding memoryless channels, for 
which quantum coding theorems hold [Hj. In particular, the quantum capacity Q 
is linijv^oo Qn/N. Equation (ITTj) by itself is a sufficient condition to prove coding 
theorems. Therefore, in the following we will use the wording forgetful channel for any 
system satisfying inequality (11), independently of the model from which memory arises. 
Now we focus on two specific, physically significant models. 

3.2. Markovian model 

The first model is a quantum channel that maps an arbitrary iV-qubit input state p 
onto 

p' = S N {p)= A ii-iNp A \ l ...i N , ik = 0,z, (12) 

ii,...,ift 

where the Kraus operators A^...^ are defined in terms of the Pauli operators cro = 1 
and a z : 

Ai...ijv = y/Ph...i N Bi 1 ...i N , Bh...i N = 0~\^ ® ■ ■ ■ (g) 0~\ N \ (13) 
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with J2{i k \Ph...i N = 1 and acting on the fc-th qubity]. 

The quantity Pi x ...i N can be interpreted as the probability that the ordered sequence 

a£\ of Pauli operators is applied to the N qubits crossing the channel. We 

define the single-qubit marginal probability p iq = J2{i h ,k^q}Ph-i N an d similarly the 

two-qubit marginal probability Pi q ,i q and assume that {pi q } = {1 — p z ,p z } for all 

q = 1, . . . , N. Under these conditions the maximum of coherent information in model 

(|T2|) is obtained for the totally unpolarized input state p unp = (l/2 N )t® N . To prove 

this statement, we construct the same iterative transformation as in ([7]) but with 
E (A0 = t (i) g . . . g t (k-i) g a (k) g t (fc+i) g . . . g t (N) instead of and notice that 

Pn — Punp is obtained starting from an input state po diagonal in the reference basis. 
Moreover it can be proven that in this case 

S[S N (^p ^)] = S[E^ Po E^\ = S( Po ). 

Since po is diagonal and £n only changes off-diagonal matrix elements, then En(po) = Po 
and S[En{^^ Po^x )) = S[£n(Po)]- We can also prove that 

S[E N ^i k) Po^ x k) )] = S[t^E N ( Po )t^} = S[E N (pa)}. 

Here S z is defined as S 2 but acts on the environment. Therefore, I c (£n, S^po^* 1 ) = 
I c (Sn,Po)- Taking again advantage of the concavity of coherent information for 
degradable channels, we finally obtain 

Ic{£N,Punp) > Ic(£n,Po)- (14) 

We can explicitly compute the quantum capacity when the joint probabilities in 
equation (1131) are described by a Markov chain [6j ITT] : 

Pii,...,i N Pi\P%i\i\ ' ' 'Pijv|ijv— 1> (1^) 

where 

Pi k \ik-i = (!- V)Pik+ V S h,ik-i- ( 16 ) 
Here p G [0, 1] measures the partial memory of the channel: it is the probability that 
the same operator (either 1 or o~ z ) is applied for two consecutive uses of the channel, 
whereas 1 — p is the probability that the two operators are uncorrelated. The limiting 
cases p = and p = 1 correspond to memoryless channels and channels with perfect 
memory, respectively. In this noise model p might depend on the time interval between 
two consecutive channel uses. If the two qubits are sent at a time interval r r c , 
where r c denotes the characteristic memory time scale for the environment, then the 
same operator is applied to both qubits (p = 1), while the opposite limit corresponds 
to the memoryless case (p = 0). 

X The Kraus operators (p~3|) define a generalized dephasing channel in the sense of equation ((4]) , with 
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The Markov chain model is forgetful, since condition (1111) is fulfilled. We first 
consider a sequence of two blocks of N + L channel uses, for which 

P = S 2 ( N +L)(p) = J2pi B iP B 1> ( 17 ) 

where the index I stands for ii, in+l+i, ^2N+l and the operators Bi are defined 
in equation ([TBI . The output state p' can be approximated by 

p' = (£ N+L r 2 (p) = Y,PiBipB\, (18) 

where the factorized probability distribution pj = Pi 1 ,...,i N Pi N+L+1 ,...,i 2N+L - Taking 
advantage of the strong convexity of trace distance [1] , we obtain 

llp'-p'lli (19) 

where the Kolmogorov distance between the probability distributions {pi} and {pi} is 
defined as 

D(pi,pi) = ^\pi-Pi\. (20) 
A i 

Using the properties of stationary Markov chains and equation (fl6l) we obtain 

D(pi,fti) < 2/i L+1 . (21) 

This implies 

\\£2(N +L )(p) - (W)) 02 (p)IIi < 2^ +1 , (22) 
from which equation ffTT]) readily followg§|. The forgetfulness of the Markov chain model 
allows us to compute the quantum capacity from the regularized coherent information 

In order to compute the quantum capacity, we consider the input state p unp and 
evaluate the coherent information I c (£n, Punp)- In this case S[SN(Pun P )] = S(p unp ) = N. 
We now take advantage of the formula (S^jv = S(W), where the density operator W 
has components W 7 ^...^,^...^ = Tr(A iu „i N pA^ t , ) [20]. Here W is diagonal and 



S(W) = - y £p {ik }log 2 p {ik} = H(X l ,...,X N ), (23) 

{ik} 

where H(Xi, ...,Xjy) is by definition the Shannon entropy of the collection of random 
variables Xi,...,X N (characterized by the joint probabilities p^...^). For a stationary 
Markov chain, we have [28] 

lim ±.H(X lt ...,X N ) = H(X 2 \X 1 )=p H(q )+p z H(q z ), 

where qo tZ = (1 — p)po )Z +p are the conditional probabilities that the channel acts on two 
subsequent qubits via the same Pauli operator, and H(q ), H(q z ) are binary Shannon 

§ It is interesting to remind the reader that the Markov chain model can also be formulated in terms 
of a structured environment [TJ1 [TJ] . 



Quantum capacity of dephasing channels with memory 




Figure 1. Plot of Qn/N as a function of fi, for the Markov chain model (fT5|) . with 
Po = 0.85. From bottom to top: N — 2,4,6,8,10 (black curves), N = 100, oo (grey 
curves). The dotted black line gives the memoryless quantum capacity. 



entropies, defined by H(q) = — q\og 2 1 — (1 — g)log 2 (l — q). Therefore, the quantum 
capacity is given by 

Q = l-p H(q )-p z H(q z ). (24) 

It is interesting to point out that Q increases for increasing degree of memory of the 
channel. In particular, for \i = we recover the capacity Q = Q\ — 1 — H(p ) of the 
memoryless dephasing channel, while for perfect memory (/i = 1) Q = 1, that is, the 
channel is asymptotically noiseless [12]. We also note that the right hand side of (1241) is 
known [11] to be a lower bound for the quantum capacity of the Markov chain dephasing 
channel. Our results prove that this bound is tight. 

In order to illustrate the convergence of Qn/N to its limiting value Q, we first 
compute the entropy exchange for the iV-qubit input state p unp . It is easy to check that 

(S e ) N = poH(qo) + p z H{q z ) + {S e ) N ^. (25) 
Using this recurrence relation we obtain 

(S e ) N = (N- l)[p H(q )+p z H(q z )] + (S e ) u (26) 
where (5 e )i = H(p ). Therefore 

Q N = N-(N- l)[p H(q )+p z H(q z )] - H(p ). (27) 

A plot of Qn/N for various iV as a function of the memory factor p, is shown in figured! 
It is clear that the convergence of Qn/N is faster when the memory factor is smaller. 
Indeed, it is easy to prove that 

tN^Q-^f (28) 

is a growing function of p, with €jy(p = 0) = and ejy(p = 1) = H(p )/N. Moreover, 
for p 1 we obtain 



e N (p) 



1 p 2 
21n2 N' 



(29) 
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3.3. Spin-boson model 

The second model of dephasing channel is defined by the system (qubits)-environment 
Hamiltonian 

H(t) = H E - l -X E F{t) + H c . (30) 

Here H E = ^ a oj a b\p a is a bosonic bath and X E = + &a) is the environment 

operator coupled to the qubits. The k-th qubit has a switchable coupling to the 
environment via its Pauli operator ay>: 

N 

F(f) = A 5>« /*(*), (31) 

k=l 

where fk{t) = 1 when the qubit is inside the channel, and /&(£) = otherwise. Finally, 

\ 2 N 

Hc = J2£rT,4 k) (32) 

a 4U; « k=l 

is a counterterm [29]. We call t p the time each carrier takes to cross the channel and 
r the time interval that separates two consecutive qubits entering the channel. The 
Hamiltonian fl30l is expressed in the interaction picture with respect to the qubits. If 
initially the system and the environment are not entangled, the state of the system at 
time t is given by the map flTJ) where 

U(t) = Te-tfi daHi ' ) . (33) 

In particular, we are interested in the final state p' = p(t = r N ), where = t p +(N — \)t 
is the transit time across the channel for the iV-qubit train. To treat this problem we 
choose the factorized basis states {|ja E )}, where - as above - {\j) = ■■■,Jn)} are 
the eigenvectors of rifc^i ■ The dynamics preserves the qubit configuration \j) and 
therefore the evolution operator (1331) is diagonal in the system indices: 

(ja E \U(t)\la' E ) = S jt (a E \U(t\j)K), (34) 

where U(t\j) = (j\U(t)\j) expresses the conditional evolution operator of the 
environment alone. Therefore 

(P% = (p);zE<«b| U{t\j)wU\t\l) \a E ). (35) 

a 

In this basis representation the environment only changes the off-diagonal elements of 
p, while populations are preserved. If the environment is initially in the pure state 
wo = |0)e(0|, then the equations (H])-© are recovered. At any rate, it is sufficient 
to consider a purification of w in an enlarged Hilbert space to write our model as a 
generalized dephasing channel (jlj). 

For a multimode environment of oscillators initially at thermal equilibrium, w = 
exp(— f3H E ), we obtain 

J2(a,\U(t\j)wU\t\l)\a E ) = 



cxp 



J IT Ld z 1 f~. 

n k=l 
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where S(w) is the power spectrum of the coupling operator X E . 

A central question is if and under which conditions a spin-boson environment 
gives a forgetful channel. Even though we cannot give a rigorous proof, we 
conjecture on physical grounds that an exponential time decay of the bath symmetrized 
autocorrelation function C{t) = 1/2 (X E (t)X E (0) + X E (0)X E (t)) is a sufficient condition 
for forgetfulness. To support this conjecture, we proof inequality ffTTj) in the particular 
case in which two single channel uses (N = 1) are separated by idle times Lt. We 
consider two qubits (M = 2 in equation (fTTl)). prepared in a generic input state p. Then 
we compute the output state p' from equation (|35|) . that is, taking into account memory 
effects, and the output p' in the memoryless limit. We obtain, for a generic monotonic 
decaying autocorrelation function, 



where the dephasing factor g is such that (p')oi = 5'(p)oi and is readily derived 
from ( 1361) by letting N = 1. In particular we consider a Lorentzian power spectrum 
S(u>) = 2t c /[1 + (cot c ) 2 }. In this case, the autocorrelation function is C(r) = e~ T / Tc and 
equation (1371) is replaced by 



Inequality fl38|) is the (fTTl) in the particular case N = 1 and M = 2 (we can set 

h = AX 2 g 2 T 2 by noting that (1 — e^ Tp ^ Tc ) 2 < 1). We conjecture that ( fTTl) also holds for 

any N and M, since the correlations between blocks of N qubits decay exponentially 

with the delay time Lt. 

A remarkable feature of model (1301) is that in the limit of perfect memory (r c — > oo) 

( f) 

there exists for any number iV of qubits a decoherence-free subspace rifj , corresponding 
to a qubit train with an equal number of |0) and |1) states. Since the dimension d of 
this subspace is such that log 2 d pa N — l/21og 2 iV at large N, then the channel is 
asymptotically noiseless, that is, Q — 1. A coding strategy naturally appears when 

blocks of N ^> 1 qubits can be sent within the memory time scale r c : if the quantum 

( f) 

information is encoded in the decoherence-protected subspace in such a way that 
the input state p is maximally mixed within this subspace, then a lower bound for 
the coherent information can be estimated as I c (£n, p)/N pa \og 2 [dim(H^)]/N pa 
1 — log 2 N/(2N). The memoryless dephasing channel instead is recovered in the limit 
t c —>■ and in this case the coherent information is maximized by the totally unpolarized 
input states p unp and the channel capacity Q — Qi — 1 — H(p ), where p = (1 + g)/2. 

Even though we could not compute the channel capacity for generic values of 
t,t p , and r c , we show in figure [2] numerical results of the coherent information I c for a 
Lorentzian power spectrum S(u) and for the input state p unp as a function of the degree 
of memory of the channel, measured by the parameter £ = t c /(t + r c ). We fix t c ,t p 
and vary r, so that the memoryless and perfect memory limits correspond to £ — > 
(r — > oo) and £ — > 1 (r — > 0). The curves in figure [2] show that memory effects enhance 
the coherent information I c /N and that I c /N grows monotonously with N . Furthermore, 
these numerical data strongly suggest that I c /N converges, for iV — > oo, to a limiting 



p'-p'h<4X 2 g 2 T 2 C(LT), 



(37) 



P ~ P'\\l < 4\ 2 g 2 T 2 {l - e -rv/r c f e -Lr/r c _ 



(38) 
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Figure 2. Plot o£I c /N as a function of £, for the bosonic bath model (|30|) : Lorentzian 
power spectrum, A = 1, r c = 1, t p = r c , maximally mixed input state. From bottom 
to top: TV = 2, 4, 6, 8, 10. The dotted line gives the memoryless quantum capacity. 

value larger than the memoryless capacity Q\. This value would provide, assuming the 
above conjectured forgetfulness for the model, a lower bound for the quantum capacity. 
Therefore, using the previously mentioned double blocking strategy, it is possible to 
increase the transmission rate if the quantum information is encoded in arbitrarily long 
blocks, separated by time intervals larger than r c . 

4. Conclusion 

In summary, we have shown that the coherent information in a dephasing channel with 
memory is maximized by separable input states, computed the quantum capacity Q 
for a Markov chain noise model and suggested a numerical lower bound for Q in the 
case of a bosonic bath where memory effects decay exponentially with time. These 
results also rely on the concept of forgetfulness, which we prove for the first model and 
strongly support on physical grounds for the second one. It would be relevant to further 
clarify the connection between the decay of environment autocorrelation functions and 
forgetfulness. It is important to point out that differently from previous works on 
quantum memory channels [6], we have carried out the limit in which the number 
of channel uses N — > oo. It would be interesting to investigate to what extent the 
results presented in this work could be applied to other physically relevant degradable 
noise models such as the amplitude damping channel [30]. Another physically relevant 
question is whether our results could be generalized to environments with algebraically 
decaying memory effects, which may model typical low-frequency noise in the solid state. 

Note: After completion of our work we became aware of a related paper [31], in 
which, in particular, the quantum capacity of a Markov chain dephasing channel is 
provided. Their derivation, not reported in that paper, is based on a method different 
from ours (S. Virmani and M. Plenio, private communication). 
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